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We outline a theory of symmetry protected topological phases of one-dimensional quantum walks. 

We assume spectral gaps around the symmetry-distinguished points +1 and -1, in which only dis¬ 
crete eigenvalues are allowed. The phase classification by integer or binary indices extends the 
classification known for translation invariant systems in terms of their band structure. However, our 
theory requires no translation invariance whatsoever, and the indices we define in this general setting 
are invariant under arbitrary symmetric local perturbations, even those that cannot be continuously 
contracted to the identity. More precisely we define two indices for every walk, characterizing the 
behavior far to the right and far to the left, respectively. Their sum is a lower bound on the number 
of eigenstates at +1 and -1. For a translation invariant system the indices add up to zero, so one 
of them already characterizes the phase. By joining two bulk phases with different indices we get 
a walk in which the right and left indices no longer cancel, so the theory predicts bound states at 
+ 1 or -1. This is a rigorous statement of bulk-edge correspondence. The results also apply to the 
Hamiltonian case with a single gap at zero. 

PACS numbers: 05.60.Gg 03.65.Db 03.65.Vf 


INTRODUCTION 


Topological phases play an important role in the classification of quantum matter, e.g. in the distinction between 
topological and ordinary insulators in lattice systems G3I3II]. Such phases are intrinsic properties of physical systems 
exhibiting (a set of) symmetries (for a review, see [29] and references therein). A key feature of topological insulators 
is, that in contrast to ordinary phase transitions covered by the so-called Landau theory [44], the symmetries of the 
system remain unbroken during transitions between distinct topological phases. 

Over the past years, such phenomena gained a lot of attraction both in theory and experiments due to their robust¬ 
ness against local perturbations which opens a wide range of applications from spintronics to topological quantum 
computation m m ns m- 

A key intuition about topological phases is a principle known as the bulk-edge correspondence. Loosely speaking 
it states that if two systems in distinct phases are joined, a bound state should emerge near the interface. Moreover, 
this should be true irrespective of how the systems are joined. This additional stability distinguishes topologically 
protected bound states from a much more commonplace phenomenon: Whenever a system is perturbed locally, for 
example by introducing an impurity or defect, we know that bound states might appear m ■ Hence even if two 
systems in the same topological phase are joined, bound states will typically arise, but these will not be topologically 
stable and can be eliminated by engineering the crossover differently. 

Classification of topological phases is well understood for translation invariant bulk systems. In that case the 
topological properties can be stated in terms of the winding properties of the energy bands over the quasi-momentum 
space, technically expressed by the K-theory of vector bundles |36| . However, this structure is lost as soon as a 
phase boundary is introduced. Hence, as remarked in [35], any theory of bulk-edge correspondence has to deal with 
non-translation invariant systems from the outset. It also has to deal with a vastly larger set of perturbations: In 
the translation invariant case, finitely many parameters suffice to characterize all systems with given maximal step 
size, but without translation invariance we obviously have to deal with infinitely many perturbation parameters. In 
addition to small global perturbations, we can have local perturbations, which cannot be achieved in many small 
steps, and are hence not amenable to standard deformation arguments [35] Sect. IV.E], It is not a priori clear that a 
meaningful classification stable under all such perturbations exists. But, as we show in this paper, the combination 
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of a gap condition with local symmetry sufficiently tames the possible behaviours at infinity to allow such a phase 
classification. It even turns out that this classification in the general case is the same as for the bulk: Indeed, every 
class in our classification can be realized by joining two translation invariant systems, and the pair of their bulk indices 
characterizes the class. This is the strongest statement of bulk-edge correspondence one could hope for. 

The aim of our paper is to describe the basics of a rigorous theory of bulk-edge correspondence in this sense, and to 
sketch the key ideas of a proof, of which the details will appear elsewhere m- The setting chosen is that of discrete¬ 
time evolutions of single-particles with internal degrees of freedom, so-called “quantum walks” 0 0 El m m HE 
As such, quantum walks have recently attracted much attention as a computational resource 0123 H3 HE \M [521- 
in particular quantum walks have been shown to exhibit a rich variety of quantum effects such as Landau-Zener 
tunneling EE the Klein paradox g3], Bloch oscillations |24j . By taking into account on-site interactions between 
two particles performing a quantum walk, the formation of molecules has been established 0 . Local defects usually 
generate bound states HE Eigenstates can also be generated by i.i.d. random 0EH] and quasi-periodic coin 

choice leading to Anderson localization or by suitably modified rules to give infinitely degenerate subspaces without 
transport (Konno localization |301 JT] ). 

Quantum walks have been experimentally realized in such diverse physical systems as neutral atoms in optical 
lattices |31j , trapped ions [55 , U’> , wave guide lattices nom] and light pulses in optical fibres |54l 5Sj as well as single 
photons in free space [T4j . 

The quantum walks considered in this manuscript obey a suitable subset of discrete symmetries { 77 , r, 7 }, commonly 
referred to as “particle-hole”, “time-reversal” and “chiral” symmetry. It is well known from the theory of bulk phases 
that in one dimension such symmetries must be imposed on the system as well as all its perturbations [22, 5Bj . In 
accordance with this theory we select the subset of the symmetries of the tenfold way [7] for which non-trivial phases 
are expected. In contrast to the ’periodic table’ in |[36] we define and compute the invariants for the various symmetry 
types in an elementary group theoretical way. 

Remarkably, our results also apply to the continuous time, i.e., (static) Hamiltonian case. However, we chose chiefly 
the discrete time setting, which applies to periodically driven “Floquet” systems, because, on one hand, we plan such 
an experimental realization in an optical lattice with neutral atoms [6], and on the other the discrete time case has 
some additional intricacies: First, there are two symmetry-related spectral gaps, at ±1 for the unitary transition 
operator, instead of the single gap at 0 in the Hamiltonian case, and secondly there are local perturbations which 
cannot be contracted to the identity while preserving the symmetry. In order to achieve a classification, which is 
invariant also under non-contractible local perturbations, new methods beyond the usual deformation arguments are 
required. 

An important input for our work were the publications of Kitagawa et al. US®] and subsequent work by Asboth 
mm- These papers give examples which realize many possible bulk symmetry classes. Supporting evidence for 
the bulk-edge principle is given by numerically finding the dynamical signature of bound states in joined systems and 
observing that they appear exactly in situations in which the bulk-edge principle predicts them UDIEM®]. 

In this note we prove this observation in full generality. Naturally, with the benefit of a precise mathematical 
formulation, we can also sharpen some of the claims and statements in the aforementioned papers. We caution the 
reader, however, that our proof so far covers only 4 of the 5 symmetry types considered in EE and is confined to the 
ID lattice case. 


BASIC FEATURES 


Let us outline some basic features of our theory. First of all we consider infinite one-dimensional systems only, but 
make no assumption of translation invariance whatsoever. Lattice sites are labelled by integers x, and the variables at 
that site are given by a finite dimensional Hilbert space 7~L X . In the translationally invariant systems all the T~L X must 
be the same and equal to the “coin space” of the walk, typically two-dimensional. The walk is a unitary operator W 
on the total Hilbert space H = ©xez^ x ’ with a finite interaction length L. That is, matrix elements between sites 
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FIG. 1. (a) schematic spectrum of unitary quantum walk with band spectrum (green) and discrete spectrum (red dots) and 
essential gaps at ±1. The arrow symbolizes the action of the symmetry operators 77,7 on quasi-energies. We are interested 
in the eigenvalues on the symmetry axis indicated, around which we assume an essential gap. (b) spectrum of self-adjoint 
Hamiltonian marked analogously. 


further apart than some finite L vanish. The discrete symmetries { 77 , r, 7 } act on each T-L x separately and we consider 
only symmetry types forcing the spectrum of the walk to be symmetric with respect to the real axis (see Fig. |Tj). The 
two real points ±1 in the spectrum of the walk IF are therefore special. These points are assumed to be in the gap of 
the bulk systems. Also it is here that the bulk-edge principle predicts protected bound states. We cover both cases 
by assuming essential gaps at ± 1 , which means that in a small interval around ±1, W has at most finitely many 
eigenvalues (counting degeneracy). This property is automatically preserved under arbitrary local perturbations [60]. 

To any walk in this setting (possibly defined only on a half-infinite lattice) we associate an index si(IF) which takes 
either integer values or the values 0,1 depending on the symmetry type considered. It is defined in terms of the 
actions of the symmetry operations in the eigenspaces at ±1 (see below). In particular, |si(IF)| is a lower bound on 
the total dimension of the two eigenspaces at ±1. For walks on the doubly infinite lattice, this quantity splits into 
two contributions 

si(IF) = si (IF) + si (IF). (1) 

The term Hi(IF) (resp. si (IF)) only depends on the behavior of IF far to the right (resp. the left). Hence these are 
independent of any local perturbations, and the pair (si (IF), si (IF)) is the basic invariant of our theory. 


A translation invariant system cannot have any isolated eigenvalues of finite degeneracy, so si(IF) = 0, and si (IF) 

already contains all classifying information. Now consider the typical scenario of bulk-edge correspondence: We take 

two translation invariant walks IFi, Wo, and a crossover between the two, i.e., a walk IF which coincides with IFi far 

-4 ■ «- -4 

to the left and with W 2 far to the right. Then we get at least | si(IF) | = |si(IF 2 ) + si(lFi)| = |si(W 2 ) — si (IFi)| 
independent eigenvectors at ±1. In particular, we get at least one if the classifying invariants are different, i.e., 
si (IF 2 ) 7 ^ si (IFi). All these statements also hold for the Hamiltonian case. 


SYMMETRY TYPES 


In every model we fix a symmetry type by specifying which of the symmetries 77 , r, 7 are present, and also the signs 
in 77 2 = ±1 and r 2 = ±1. The particle hole symmetry 77 and the time reversal symmetry r are always antiunitary 
E 2 E 2 , and the chiral symmetry 7 is unitary. When two such symmetries are present, their product provides the 
third, and we choose the phase convention making them commute. The three signs in Table |T| then multiply to 1. 
A Hamiltonian H is said to satisfy the symmetry if rjH = — Hr 7 , yTT = —H 7 , and tH = Hr. The corresponding 
conditions for walks can be read off the conditions which the time evolution operator IF = e~ ztH satisfies in this case, 
i.e., 77 IF = IF 77 , 7 IF = IF* 7 , and rlF = IF*t, if the respective symmetry is present. The symmetries are assumed to 
act in each cell 'H x separately. We also assume that the symmetry operators in each cell are balanced meaning that 
there is some unitary acting only in the cell, which satisfies the symmetry and has no eigenvalues ± 1 . 

Given the symmetry operators we now consider quantum walks IF which (1) satisfy the symmetry and (2) are 
essentially gapped. These assumptions carry over in an obvious manner to the Hamiltonian case. 
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TABLE I. Symmetry types considered in this paper, with their generators and relations for their squares. Absence of an entry 
in the respective column means that the generator is not part of the type. The second column gives the Cartan classification [7]. 
The perhaps unusual sign of y 2 in the last row is due to the convention which makes the three symmetry operators commute. 
I (S) is the range of the symmetry index si of a finite dimensional representation of the symmetry. The last column gives an 
explicit expression for si, where d denotes the dimension of the representation. 


DEFINING THE SYMMETRY INDEX 


The main idea for the definition of the symmetry index si(W) is now suggested by perturbation theory. We want 
the index to be stable under any small deformation of W which respects the symmetry. As an example, consider the 
first symmetry type (77 only). Then, since eigenvalues 7 ^ ±1 occur in pairs, the dimension of the eigenspace at +1 can 
change, but only by multiples of 2. The same holds for the eigenspace at —1 and hence for the total dimension. We 
thus take si(W) = 0,1 as the parity of the combined dimensions of the eigenspaces at ± 1 . 

We can generalize this to all the other symmetry types by observing that the change of the eigenspaces can only be 
by a subspace, on which the symmetries act in a balanced way as defined in the previous section. Concretely, for the 
types 2,3,4 in the table, consider an eigenvector W<f> = A cj> with 3rnA 7 ^ 0 then 7 <f> is an eigenvector with eigenvalue 
A. Therefore, on the subspace spanned by </>, 7 (f>, rjq !>, T(f>, 7 acts like a swap operator, which has zero trace. This is the 
sort of subspace by which the eigenspaces of W at ±1 may change. Therefore, if £ denotes the sum of the eigenspaces 
of W at ±1 we can take si(W) = trg 7 with the trace evaluated just over £. This ensures that si(W) is invariant under 
all small deformations satisfying the symmetry. The technical details of this argument and the appropriate version 
for the 5. symmetry type will be given in p 8 | . 


THE SPLITTING CONSTRUCTION 


The second main ingredient of our theory is a way to split a walk into a right and a left block, acting only 
on cells x < Xq and x > Xq, respectively. This is done by a unitary operator V, which differs from the identity 
only on finitely many cells around Xq, i.e., is a local perturbation. The perturbed walk is then to be of the form 
W' = VW = Wl © Wr, and must satisfy the symmetry. It is crucial that we can choose V to be a gentle perturbation , 
i.e., one that is continuously connected to the identity such that the walks connecting W and W' all satisfy the 
symmetry. We then define si (W) = si(W#) and si(W) = si (My-.). By invariance of si(W) under small perturbations, 
si(W) = si(VW) = si(Wi) + si(WA), which is ([!]). The core of the decoupling construction stems from E3 Lemma 4], 
and details will be given in |T53. We have to exclude walks of the last symmetry type in Table [I] at this point, because 
we could not establish in full generality the existence of a decoupling which satisfies the symmetry. 


CUT-INDEPENDENCE 


The topological stability of the symmetry index hinges on the observation that si(W) (and analogously si (W)) is 
independent of the cut position. To see this, consider two cut positions xq,X\ with decoupling unitaries Vo ; V) where 
x\ — Xq is sufficiently large such that the regions where these operators differ from the identity do not overlap. We 
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FIG. 2. Sketch of the decouplings needed to prove the independence of the index si (W) of the cut position. It is assumed that 
x 0 -C Xi such that the nontrivial supports of Vo and Vf do not overlap. 


allow Vq to be an arbitrary decoupling, but choose V\ to be gentle, as we may. Then V % W — Wl,i 0 Wr,i-, for * = 1,2. 
Moreover, the piece Wr$ is further decoupled by V), so that V\Wr.o = Wc 0 Wr ,i (see Fig. [2]) for some symmetric 
Wc on the cells xo, ■ ■■ ,x\ — 1. Because V\ is gentle (i.e., contractible to the identity), we have 

si(W fl , 0 ) = si(Wc) + si(W fl ,i)- ( 2 ) 

We will presently argue that si(Wc) = 0. Hence the candidates for si(W), computed for different cut positions 
coincide, and we can take the cut position x\ as large as we wish. 

It remains to consider si(Wc) for an arbitrary unitary Wc satisfying the symmetry on a finite collection of cells, 
i.e., on the space He = (BxLx 0 ^ x ■ To this end we split this space into the sum H es of the eigenspaces of Wc at 
±1 and the orthogonal complement H^ s . The symmetry operators are block diagonal in this decomposition and, e.g., 
7 c = 7es © 7i- On Hf~ s the symmetries act in a balanced way, because Wc acts as a gapped unitary satisfying the 
symmetry. Hence try^ = 0. Similarly, because on each cell, and hence on 'He the symmetry operators are balanced, 
tr 7 c = 0. Since tryc = try es + try", we get si(Rc') = tr 7 es = 0. The argument with just particle-hole symmetry is 
analogous. 


IDENTIFYING TOPOLOGICALLY PROTECTED EDGE STATES 


In the existing literature (see e.g. [12, 00]) topologically protected edge states are detected indirectly by starting 
a joined-up walk at the boundary and letting it run for a few steps. Then one looks for some probability getting 
trapped near the boundary. This can be followed nicely in a Mathematica demonstration written by Kitagawa. 

In order to rigorously pinpoint the eigenvalues we used a general method for reducing the eigenvalue problem for 
a locally perturbed infinite system to that finite subsystem on which the perturbation happens. The influence of the 
infinite system is then taken into account by the matrix valued Schur function of the unperturbed system (which 
is related to the Green function). Some background on Schur functions can be found in [T51 ITQI ESI and references 
therein. For any unitary W on H and any finite dimensional subspace Ho C H the Schur function f(z) is an analytic 
function on the open unit disc with values in the operators on Ho , satisfying ||/(z)|| < 1. In an essential gap / is also 
analytic and unitary on the unit circle. The eigenspace of W for eigenvalue z, or rather the eigen-subspace in the span 
of all W n Ho is then canonically isomorphic to the fixed point space of the unitary zf(z). This is ideally suited for 
perturbations on finite dimensional subspaces, because if V is such a unitary differing from 1 only on Ho, the Schur 
function of VW is simply fvw{z) = f(z)V* [T9] . Moreover, if Ho is chosen consistent with a discrete symmetry, this 
whole construction is also consistent with the symmetry, so one can immediately read off the eigenspaces and the 
relevant representations of the symmetry, and hence obtain si±(V"W). 

The Schur function of a quantum walk is readily calculated using the CGMV approach to quantum walks mnzi 
mm, which proceeds by reducing the unitary evolution to canonical form, given by the so-called CMV matrices 
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nangnziEBi. Moreover, by using for V the splitting around the origin introduced above one obtains a half-space 
walk in CGMV form, which coincides with the given, translation invariant one to the the right of the origin. One 
thus only needs the Schur function relative to the subspace on which V differs from the identity, i.e. at finitely many 
cells around x = 0. From this the exact phase diagram (e.g. Fig. [3] (left)) follows. 


LOCAL VS. GENTLE PERTURBATIONS 


As Fig. [T] shows one difference between the unitary case considered here and the the Hamiltonian case is that we 
have two gaps. It is suggestive to try to associate invariants with each of these gaps El. since the basic perturbation 
argument showing si(IT) to be constant under symmetry preserving deformations applies to both gaps separately. So 
let us introduce si+(lT) as the symmetry index of the symmetries acting in the eigenspace of W for eigenvalue +1, 
and similarly si_(W) for the — 1-eigenspace. Combined with a splitting VW = Wr © Wr we get the four numbers 
in the upper left of Table [IT] Summing the two columns of this 2 x 2 -matrix yields the invariants we have already 
defined. They are insensitive to any local modification of the walk, and in particular independent of how the cut is 
made. They are also invariant under arbitrary norm-small perturbations CHI, although we have not shown that here. 
Indeed, the most straightforward attempt to prove that actually fails, i.e., it is not possible [18] to find for every 
continuous family of walks t, W t gentle decouplings V t which also depend continuously on t. On the other hand, 
if the decoupling is chosen to be gentle, we have that Wl © Wr is a continuous deformation of W, and so each row 
adds up to si±(W). The gentleness of the decoupling is a crucial condition here: There are also non-gentle, but local 
decouplings, and for these the sum rule does not hold. Finally, the individual numbers si±(Wi, j .R) have no useful 
invariance properties at all: they are neither invariant under local nor gentle perturbations, and generally depend on 
the position of the decoupling cut. An example of the indices si + (IT) and si_(IT) depending on the type of the cut 
is given in Figure [4] 


si + (Wi,) 

si+(ITfl) 

si+(IT) 

bL. (Wl) 

si-(ITfl) 

si (IT) 

si (IT) 

si (IT) 

si(IT) 


TABLE II. Table of four potential invariants from the splitting construction. The sum in each row gives the marginal shown, 
provided the decoupling is chosen to be gentle. The column sums are always the invariants discussed in this paper. Individual 
entries have no invariance properties, but the row sums are invariant under all gentle perturbations. 


To summarize, we have two invariants (si (IT), si (IT)) for all perturbations, but if we only discuss gentle pertur¬ 
bations, there is a third invariant, namely si+(IT), i.e., the sum of the first row in Table |ll| The second row would 
give the same information, since si_(IT) = si (IT) + si (IT) — si+(W). 


EXAMPLE: THE SPLIT-STEP WALK REVISITED 


Introduced in m, this example is also treated in pm HU ESI ES], and many other papers. It has symmetry type 3, 
see Table ^ The coin space is two-dimensional, with rj the complex conjugation in position space. Its name derives 
from the use of two separate shift operations, S)-, the right shift of the spin-up vectors and the left shift of the 
spin-down vectors. Choosing the local basis in each cell 7 ~L X appropriately (sometimes called “choosing a different 
time-frame” ED the split-step walk takes the form 


IT = BS^AS^B, (3) 

where A = (J),. A x and B = B x are unitary operators acting sitewise, and each A x , B x is an admissible unitary 
operator acting on 7 i x . In the basis chosen the chiral symmetry takes the standard form 7 = <Ji and 77 , if 
applicable, is given by complex conjugation. In the form of © it is straightforward to see that admissibility of the 
A x , B x suffices for the admissibility of IT. To see this, note that admissibility for 77 depends solely on the reality of 
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the entries of W and hence follows by choosing A x = R(6 2 ) and A x = R{9\/2) to be real rotations since Sf and S 4 , 
are real anyhow. For the chiral symmetry 7 , we use the indentities 

7<Sl = s; 7 , 7-St = s tr (4) 

to prove by straightforward computation that hermiticity of W 7 follows directly from the hermiticity of B 7 and Ay. 
The phase diagram for the index si (W) is shown in Fig. pi 



FIG. 3. Interactive version under [T|. 

Left: Phase diagram showing the range of si (IV) for translation invariant split-step walks with angle parameters 61 , 62 . The 
white lines represent gap closures, hence systems not covered in our setting. Right: Chirally off-diagonal matrix element of 
walk parameterized by quasi-momentum. The winding number of this curve around the origin is the topological index for 
translation invariant systems. This bulk index coincides with the symmetry index si (IV), which is shown explicitly for the 
three points A, B, C. 

This approach to split-step walks allows one to easily generate examples with phase diagrams of arbitrary dimension. 
The symmetric form suggests a straightforward extension to a split-step walk with three-dimensional phase diagrams 
by considering 

W = CSfBSfAS+BS+C. (5) 

The admissibility of this walk follows on the same lines as for ([3]). To generate an example of symmetry type 2 but 
explicitly violates particle-hole symmetry r/ the rotation matrices C were garnished by C 1 —>■ C (_\ ) /y/ 2. The 

chirally off-diagonal matrix element as well as the eigenfunctions are plotted in Figure [4] The specific angles of the 
rotations A, B , C lead to a winding number si (IV) = +1. 


SUMMARY AND OUTLOOK 


We have established a theory of topological invariants for one-dimensional quantum walks with discrete symmetries. 
The known classification of translation invariant phases is reproduced, and in addition we get a full classification of 
non-translation invariant systems stable under arbitrary symmetric local perturbations, thus providing a rigorous 
theory of bulk-edge correspondence in this setting. The technical details will be presented in a future publication [li§| . 

We believe that some variant of this theory will also apply to the fifth symmetry type in Table [TJ which so far is 
not covered. Moreover, we envisage extensions to higher spatial dimensions and the periodic table 3b of symmetric 
systems and edge modes, at least when the translation symmetry parallel to the interface between two phases is not 
broken. 
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FIG. 4. Left and Center: Real parts of eigenfunctions of the decoupled walk W. The red color encodes chirality —1, whereas 
blue indicates +1. A solid line is used when the corresponding eigenvalue of W is 1, where the dashed line indicates —1. The 
decoupling matrix at x = 0 is chosen as both, Aq and Bo being the reflection matrix g\ (left) contrasted to the decoupling being 
realized by both, Ao and Bo being the gentle decoupling iao. Clearly visible in both cases is the localization on either the left 
or the right side of the cut position. The insets reflect Table [ll] and clearly show that whereas both, si (IT) as well as "st(IT) 
are invariant under the decoupling chosen, the homotopy invariants si_(IT) and si + (lT) are not invariant under non-gentle 
perturbations, thus confirming the abstract results above. 

Right: Typical winding of chirally off-diagonal matrix element of leading to si (IT) = +1. 
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